Constraining nonextensive statistics with plasma oscillation data 
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We discuss experimental constraints on the free parameter of the nonextensive kinetic theory from 
measurements of the thermal dispersion relation in a collisionless plasma. For electrostatic plane- 
wave propagation, we show through a statistical analysis that a good agreement between theory 
and experiment is possible if the allowed values of the g-parameter are restricted by q — 0.77 ± 0.03 
at 95% confidence level (or equivalently, 2 — q = 1.23, in the largely adopted convention for the 
entropy index q). Such a result rules out (by a large statistical margin) the standard Bohm-Gross 
dispersion relation which is derived assuming that the stationary Maxwellian distribution (q — 1) is 
the unperturbed solution. 
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It is usually assumed that the particle velocity distribu- 
tion in a collisionless thermal plasma eventually relaxes 
to the standard Maxwellian curve 0,0. However, the 
possibility of a non-Gaussian behavior was established 
since the earlier experiments in plasma physics Q . More 
recently, Langmuir probe measurements have also sug- 
gested that the isotropic component of the electron dis- 
tribution resulting from inverse bremsstrahlung absorp- 
tion in plasma is well described by a non Maxwellian 
distribution 0. In principle, such results may be ex- 
plained by a stationary distribution that emerges nat- 
urally from the nonextensive statistical formalism pro- 
posed by Tsallis Such an approach generalizes the 
standard Boltzmann-Gibbs formalism through a new an- 
alytic form for the entropy S q = fce(l — YliPi)/(l ~ 1)' 
where fee is the standard Boltzmann constant, pi is the 
probability of the ith microstate, and q is a parameter 
quantifying the degree of nonextensivity. This expres- 
sion has been introduced in order to extend the applica- 
bility of statistical mechanics to system with long range 
interactions and has the standard Gibbs-Jaynes-Shannon 
entropy as a particular limiting case (q = 1). The asso- 
ciated nonextensive kinetic theory has also been consid- 
ered in many different physical scenarios, ranging from 
astrophysics to plasma physics (see @ for a regularly 
updated bibliography). For example, it has been suc- 
cessfully applied to two-dimensional Euler and drift tur- 
bulence in a pure-electron plasma column Q , Levy- type 
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anomalous diffusion 8], anomalous relaxation through 
electron-phonon interaction 0, ferrofluid-like systems 
[To| , plasma oscillations the solar neutrino prob- 

lem [l^, astrophysical systems [l^, among others. This 
new formalism has also shown to be endowed with sev- 
eral interesting mathematical properties [T^j. with the 
main theorems of the standard statistics admitting suit- 
able generalizations More recently, some authors 
have argued that the value of the nonextensive param- 
eter is predictable for some simple physical systems, as 
for instance, in the case of logistic maps. Their ana- 
lyzes take into account the relation between sensitivity to 
initial conditions and relaxation with the q-generalized 
Lyapunov coefficient usually computed in terms of the 
map parameters [Tiij ]. Theoretically, the situation is not 
so neat for long-range Hamiltonian systems as it is for 
maps. The present belief is that systems with short 
range interactions obeys the Boltzmann-Gibbs statistics 
and exponential relaxation dictates the approach to equi- 
librium. However, for more complex systems (including 
the case of plasmas) , it is not possible to determine apri- 
ori how the relaxation to a stationary state takes place. 
In some cases, due to a formation of correlated clusters 
in N-Body systems, a preequilibrium stage may happens 
much earlier than the final Boltzmann-Gibbs equilibrium 
in such a way that the overall effect can be described 
by a power-law Tsallis distribution ^j- The interesting 
feature of this power-law function is that many models 
are analytically tractable so that a detailed comparison 
with the Boltzmann-Gibbs approach is immediate. In 
this letter, we discuss new constraints on the g-parameter 
associated to measurements of the dispersion relations 
for electrostatic plasma oscillations. Through a statis- 
tical analysis we estimate the value of the nonextensive 
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parameter from the existing experimental data. More 
precisely, by analyzing the data set for thermal disper- 
sion relation taken from Van Hoven |l8j we argue that 
a reasonable agreement between theory and experiment 
is possible for values of the nonextensive parameter re- 
stricted by q = 0.77 ±0.03 at 95% confidence level. Some 
years ago, it was shown how Maxwell's derivation ^| 
to the stationary velocity distribution function could be 
extended to the nonextensive domain |2(J, |2lJ. In the 
one-dimensional case, the new distribution reads 



fo(v x ) = A q 



l-(3-l) 



2k B T 



(1) 



The above expression has been determined from two 
simple requirements: (i) the isotropy of the veloc- 
ity space, and (ii) a nonextensive generalization of 
the Maxwell factorizability condition, that is, F(v) ^ 
f( v x)f(v y )f(v z ) 20]. The quantity A q in the above equa- 
tion denotes the ^-normalization constant and can be 
written, respectively, for the intervals < q < 1 and 
Q > 1. by 



A, 



0<q<l 



nT(a) / m(l — q) 
r(a-i)V 2nk B T 



and 



A q>x = 
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where a = l/(q — 1) is a dimensionless number, n is the 
particle number density, m is the mass of the particle 
and T is the temperature. For q > 1, the distribution 
exhibits a thermal cutoff on the maximum value allowed 
for the velocity of the particles. However, we notice that 
the above velocity distribution is parameterized with the 
nonextensive parameter shifted by q — > 2 — q. This means 
that our thermal cutoff at q > 1 is equivalent to the cut- 
off condition, q < 1, usually adopted by other authors. 
As one may check, in the extensive limit, q = 1, the 
above distribution function reduce to the standard Gaus- 
sian form 0, 0, |2(j • Another interesting result leading 
to a coherent nonextensive kinetic theory is related to 
the transport equation. In this context, the generalized 
Boltzmann's equation reads fgJJ 



dt 



Or 



df 
dv 



C q (f), 



(4) 



where the left-hand-side is the total time derivative of / 
while C q is the g-nonextensive collisional integral term 
responsible by changes in the distribution function /. In 
this approach, the stationary velocity g-distribution was 
obtained from a slight generalization of the kinetic Boltz- 
mann i?-theorem. The nonextensive ingredients follow 
simply by modifying the molecular chaos hypothesis and 
by generalizing the local entropy expression according to 
Tsallis' argument. One important related result is that 



the proof of the ff g -theorem is possible (in accordance 
to the second law) only if the nonextensive parameter is 
positive definite, thereby fixing the lower bound assumed 
in equation (2). On the other hand, we know that high 
frequency vibrations in a collisionless electronic plasma 
may kinetically be described in a highly simplified man- 
ner, where electron-electron and electron-ion collisions 
are unimportant, in such a way that the collisional inte- 
gral term in the Boltzmann equation may be neglected 
[l], |22| • In the linear approach, the distribution function 
of electrons is perturbed in first order, while the distri- 
bution function of ions can be considered as an invari- 
able quantity. If fo(v) corresponds to the unperturbed 
homogeneous and time-independent stationary distribu- 
tion, the resulting particle distribution function may be 
approximated as 



/ = /o(«) + /i(r,v,i), /i«/ , 



(5) 



where f\ is the corresponding perturbation in the distri- 
bution function. The dynamical behavior of the plasma 
can be described by a combination of the linearized 
Vlasov and Poisson equations. One obtains [H.l22]| 



dfi . df! e df 

~5T + v -^~ = V( ? 1 -"a — ' 

at or m ov 
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where <f>\ (r) is the first order correction in the electro- 
static potential and e denotes the electronic charge. The 
set of coupled equations and J7| may be worked out 
through the simplified derivation for electrostatic waves 
(longitudinal plasma waves), where the dispersion rela- 
tion can be calculated either by taking the constraint of 
null permitivitty or by using the mathematical tech- 
niques of integral transform originally developed by Lan- 
dau [2^ ■ In this case, the solutions of the above equations 
can be written as f\ (r, v, t) — f\ (v) exp i(k ■ r — tut) and 
4>i(r,t) = 0expi(k ■ r — cot) provided that f\ and 4> sat- 
isfy the relations 
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and 



k z (j) = Aire / /i(v)d J w 
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Finally, combining these expressions we find that the dis- 
persion relation between to and k is given by 



1 



Aire 1 
k 2 m 



k • dfo /dv 3 



<Fv = 0. 



to 



(10) 



Now, by considering the limit case of small wave number 
k <C Ud (from now on the subindex D stands for Debye 
quantities) we expand the integrand of (|10fl in power of 
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FIG. 1: Thermal dispersion relations for nonextensive velocity 
distribution. The selected values of q is showed in the picture. 
Note that nonextensive distribution with q < 1 is strongly 
suggested by these data. 



FIG. 2: x 2 — 1 plane provided by the data of Van Hoven 
[14]. The horizontal line indicates the 95% confidence limit 
for the nonextensive parameter. From this analysis we find 
q = 0.77 ±0.03 (95% c.L). 



k and substitute the nonextensive stationary distribution 
given by Q into the dispersion relation, one obtains (see 
reference [ll| for more details) 

W 2 ^l + i{\ D k) 2 ^—^, (11) 

where W — u>/u>o, lu q = ( 47r ^ e ) 1 / 2 is the natural oscilla- 
tion plasma frequency and = (k B T '/ ' Aime 2 ) 1 / 2 is the 
electronic Debye-Hiickel radius. As one may check, for 
the extensive limit 5=1, Eq. (|11|) reads 

W 2 = l + 3(X D k) 2 , (12) 

which is the standard Bohm-Gross relation pj. In or- 
der to check quantitatively the validity of this new ap- 
proach, as given by (|11|) . we consider the experimen- 
tal data set taken from Van Hoven [l8|. These data, 
originally composed by 40 experimental points, was dis- 
tributed into 6 bins by using the standard steps of data 
analysis (see, for instance, [23 )■ To find the confidence 
limits (c.l.) we use a \ 2 minimization ne glec ting the for- 
bidden region by the 7i-theorem (q < 0) [21|. i.e., 

2 (h , x ^ [W(q) - W m f 
x{k\ D ,q) = 2_^ j 1 ( 13 ) 

i=l « 

where W(q) is given by Eqs. (|llfl and (|12fl and W i is 

the experimental values of the ratio lv/uo with errors <7; 
of the i th bin in the sample. 



In Fig. 1 we show the binned data of the ratio W as 
a function of dimensionless parameter fcArj for some se- 
lected values of q. Thick solid line stands for the extensive 
Bohm-Gross result based on the Maxwellian distribution 
(q = 1). As can be seen, for values of q > 1 the curves 
increase less rapidly than in the Maxwellian case and 
clearly depart from the experimental points. In the com- 
plementary range, < q < 1, the curves increase more 
rapidly than in the Maxwellian case, showing a better 
agreement with the data. In Fig. 2 we display the result 
of our statistical analysis in the x 2 ~<Z plane. The horizon- 
tal dashed line corresponds to 95% confidence level (2er). 
For these data the peak of the likelihood is located at 
q = 0.77 with the corresponding error interval of ±0.03. 
Such a value suggests a power law behavior without ther- 
mal cutoff for the allowed values of velocities. These re- 
sults clearly show that the new formalism discussed here 
may provide a better fit for these experimental data than 
does the standard extensive approach (q = 1). Moreover, 
the present discussion reinforces the interest for new ex- 
periments involving the dispersion relations in a collision- 
less electronic plasma in order to check more accurately 
the consistency of the results predicted by the nonexten- 
sive formalism. 

In conclusion, we stress that the lack of a reasonable 
explanation to the physical meaning of the q parame- 
ter nowadays for long-range Hamiltonian systems points 
naturally to the following strategy: instead of paying 
attention to more formal results and mathematical ex- 
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tensions, it seems more important to constrain its value 
via a number of different set of experiments (see, for in- 
stance, 0), in order to determine more clearly the reality 
of nonextensive effects. Actually, it has been shown [24J 
that the manganites (material that exhibit long-range in- 
teractions and fractal geometry) can be studied through 
the nonextensive statistic. In this regard, in contrast 
with long-range interactions of manganites, the present 
study shows that the role of long-range interactions does 
not seem to be a fundamental ingredient for a statisti- 
cal analysis of thermal dispersion relation in plasma. In 
particular, we suspect that the extended kinetic theory 
can also be linked with statistical correlations through 
the velociti es ( before and after of a collision) of particles 
(see Refs. [M EH). Finally, we have showed that the 
Bohm- Gross dispersion relation which is a direct conse- 
quence of the Maxwellian distribution (q = 1) is strongly 
deprived by the experimental data, being ruled out by 



a considerably statistical margin. Our analysis indicates 
that the best fit for these particular data set occurs for 
a nonextensive distribution without thermal cutoff with 
q ~ 0.77 in our convention for the q-exponential (see 
equation (2)). Note, however, that if the standard defi- 
nition is considered (see for instance 

[HEl) these data 
set constrains the nonextensive parameter to the value 
q = 1.23 which is equivalent to a translation q — > 2 — q. 
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